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Abstract 

The  report  addresses  a  numerical  investigation  of  the  trade-off  between  some  state  space 
and  FIR  filtering  algorithms  intended  to  provide  optimal  time  error  steering  of  a  local  crystal 
clock  employing  the  GPS  reference  time.  For  the  sake  of  low-cost  crystal  clocks,  we  employ 
four  most  simple  structures,  namely  1)  Two-state  Kalman,  2)  Three-state  Kalman,  3)  FIR  with 
the  constant  kernel  (simple  MA),  and  4)  FIR  with  the  linear  kernel  (optimally  unbiased  for  lin¬ 
ear  phase  drift).  The  optimal  control  problem  is  numerically  solved  in  a  sense  of  least  mean 
squares  for  the  GPS-based  time  error  database  and  for  different  estimators  tuned  to  have  the 
same  time  constant.  We  study  a  digital  control  loop  of  a  local  clock,  assuming  the  latter  to  be 
linear  and  inertialess  with  respect  to  the  filter  memory.  We  also  compare  the  algorithms  for  the 
mean  square  error  produced  and  for  sensitivity  to  variations  of  the  feedback  coefficient. 


INTRODUCTION 

Fast  and  accurate  locking  of  a  local  clock,  along  with  accurate  estimation  of  its  performance  (time  error), 
still  remains  a  key  problem  in  GPS-based  timekeeping.  Usually,  it  is  solved  employing  two-  or  three- 
state  Kalman  filters  via  the  recursive  algorithms;  finite  impulse  response  (FIR)  filters  with  a  constant,  lin¬ 
ear,  or  exponential  kernel;  and  infinite  impulse  response  (HR)  filters  with  a  limited  number  of  poles  (as  a 
role,  no  more  than  one-two  poles).  In  our  previous  reports  [1-4],  it  has  been  numerically  shown  that  no 
one  of  the  above  listed  designs  is  universal  for  the  crystal  and  rubidium  clocks  in  terms  of  minimum  pro¬ 
duced  errors.  The  problem  may  be  distinguished  to  lie  in  two  planes: 

1.  Uncertainty  of  a  clock  time  model,  which  is  postulated  by  the  standards  to  be  of  three  states  (time 
error,  fractional  frequency  offset,  and  linear  fractional  frequency  drift  rate)  and  in  which  all  other  un¬ 
necessary  but  possibly  non-zero  states,  including  noise,  are  removed  to  the  random  phase  deviation 
component; 

2.  Different  robustness  of  filters  for  different  frequency  offsets  and  drifts  so  that  the  simplest  FIR  struc¬ 
ture  with  a  constant  kernel  becomes  best  if  the  frequency  offset  falls  below  some  threshold  [4]  and,  if 
not,  the  problem  may  be  solved  in  different  ways. 

One  more  problem  that  might  also  be  mentioned  is  that  the  clock  noise  is  inherently  colored  and,  thus,  the 
standard  Kalman  estimator  needs  to  be  modified  (Kalman  claims  the  noise  to  be  white).  Such  a  modifica- 
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tion  may  be  done  with  additional  states.  However,  in  the  GPS-based  measurement  the  receiver  Gauss 
noise  dominates  and,  therefore,  the  Kalman  filter  may  be  used  straightforwardly. 

It  has  been  a  well-known  rule  of  thumb  that,  in  contrast  to  the  HR  structure  filter,  the  FIR  one  is  more  ro¬ 
bust  against  numerical  errors  such  as  computational,  quantization,  and  round-off  Furthermore,  since  the 
FIR  filter  utilizes  finite  measurements  of  time  errors  over  the  most  resent  time  interval,  this  filter  may 
demonstrate  robustness  against  temporary  uncertainties  of  the  clock  time  model  and  even  temporary 
losses  of  GPS  signal.  In  turn,  the  HR  structure  in  such  a  situation  suffers  from  a  divergence  phenomenon. 

In  this  report,  we  numerically  investigate  the  trade-off  between  the  state  space  Kalman  algorithms,  which 
have  an  HR  structure,  and  FIR  structure  algorithms  by  utilizing  them  in  the  clock  control  loop  to  solve  the 
filtering  and  one-step  prediction  problem.  For  the  sake  of  low-cost  crystal  clocks  with  limited  memory, 
we  restrict  the  studies  by  employing  only  four  of  the  most  simple  filters,  such  as  1)  Two-state  (2D)  Kal¬ 
man,  2)  Three-state  (3D)  Kalman,  3)  FIR  with  a  constant  kernel  (CK),  and  4)  FIR  with  a  linear  kernel 
(LK).  We  model,  simulate,  and  investigate  the  GPS-based  discrete  time  control  loop  of  a  synchronized 
time  scale,  assuming  the  clock  to  be  linear  and  inertialess  with  respect  to  the  filter  memory.  We  find  nu¬ 
merically  an  optimal  loop  gain  for  each  of  the  estimators  and  evaluate  the  mean  square  error  (MSB)  of  the 
steered  time  scale.  We  finally  bring  a  numerical  example  for  the  GPS-based  digital  steering  of  a  crystal 
clock  time  error,  assuming  the  crystal  oscillator  to  be  uncontrollable. 


DIGITAL  LOOP  WITH  1-STEP  PREDICTIVE  CONTROL  OF  A 
LOCAL  CLOCK 


Figure  1  sketches  a  GPS-based  digital  control  loop  for  a  local  clock.  The  time  error  deterministic  func¬ 
tion  of  a  clock  is  extended  to  the  Taylor  series  and  traditionally  limited  to  2-3  states.  Accordingly,  its 

vector  of  3x1  or  2x1  dimensions  reads,  respectively,  or  ,  where  the  unknown 

time  error  ,  fractional  frequency  offset  ,  and  linear  fractional  frequency  drift  rate  are  calculated 
by  the  recursive  equations  for  the  model  described  by  3 -states  and  2-states,  respectively, 

3-states:  x„  =  V ,  (a)  =  y^_^  +  ,  (b)  Z)„  =  Z)„_, ,  (c) 

(1) 

2-states:  x,  =  A ,  (a)  =  y^_^ ,  (b)  (2) 


where  A  =  t^— 1^_^  is  a  sample  time^  between  two  neighboring  measurements  taken  at  discrete  time  in¬ 
stants  t^  and  t^_^ .  The  random  constituent  of  a  clock  error  suffers  of  both  the  oscillator  white  and  colored 
Gauss  noises^  [e.g.,  white  phase  (f^),  flicker  phase  (f  ^),  white  frequency  (f^),  flicker  frequency  (f^),  and 
random  frequency  walk  (f"^)],  so  that,  generally,  the  noise  vector  of  3x1  or  2x1  dimensions  reads,  respec- 
tively,  or  where  are  relative  zero-mean  colored  sta¬ 

tionary  Gauss  noises  with  known  covariances  {CxnCxm)'  {CynCyn^^  {CDnCom)  or  corresponding 

power  spectral  densities  (psd’s).  The  noises  with  different  psd  slopes  are  usually  assumed  to  be  uncorre¬ 
lated  having  zero  cross-covariances.  The  deterministic  GPS  second  (IPPS)  is  supposed  to  be  reference 

^  A  sample  time  is  a  matter  for  optimization  and  for  different  clocks  ranges  from  tens  seconds  to  1000  sec  [5]. 

^  Here  f  ,  where  /  is  integer  and  f  is  a  Fourier  frequency,  means  the  slope  of  the  certain  part  of  the  power  spectral 
density  (psd)  of  the  oscillator  phase. 


250 


35^^  Annual  Precise  Time  and  Time  Interval  (PTTI)  Meeting 


and,  thus,  error-free,  =  0 .  However,  at  the  receiver,  the  pulses  (IPPS)  are  contaminated  by  large 
zero-mean  additive  white  Gauss  noise  with  known  (measured)  intensity.  Accordingly,  the  reference 
time  error  becomes  a  pure  noise.  The  time  interval  meter  provides  measurements  of  the  clock  time  error 
employing  the  GPS  reference  signal,  and  this  error  becomes  observable  via  the  noisy  vector  . 

The  estimator  first  solves  a  stochastic  filtering  problem  to  provide  an  accurate  estimate  =  ^^i.^n\n-k) 
of  the  clock  states  (1)  or  (2)  via  .  In  this  estimate,  the  coefficient  k  is  equal  to  unity  for  the  state  space 
model,  and  it  ranges  from  1  to  A  for  the  FIR  structure.  Respectively,  an  operator  of  the  estimation 
differs.  In  the  holdover  case  when  the  GPS  timing  signals  are  not  temporary  available,  the  estimator 
forms  a  recursive  1-step  autonomous  prediction  =  3^+i(^„l„_^+i)  to  make  it  possible  for  the  algo¬ 
rithm  to  be  robust  and  for  the  clock  to  be  locked  utilizing  the  recursive  noise-free  model.  Should  the  pre¬ 
dictor  be  utilized  in  the  algorithm,  the  estimation  problem  is  crucial  for  the  clock  error  steering,  and  we 

therefore  discuss  its  solution  in  detail  in  the  following  sections.  Finally,  the  predicted  value  is 
scaled  with  an  unknown  negative  loop  gain  —  and  then  mixed  with  the  clock  1-step  ahead  actual 
value  .  The  optimal  control  law  is  then  a  linear  feedback  law  that  specifies  X^^^  to  be 

X^^^  =  ,  where  the  feedback  coefficients  matrix  for  3-states  and  2-states  is  described  by,  re¬ 

spectively. 


k 

0 

0  ■ 

xn 

0 

k„„ 

0 

and  k„  = 

kxn 

0  ■ 

yn 

n 

0 

k,,„ 

0 

0 

^Dn_ 

yn 

in  which  the  components  ,  k^^  ,  and  kp,^  are  intended  to  govern  each  of  the  clock  states. 


To  obtain  a  time-invariant  value  of  the  feedback  gain  (3)  over  the  observable  database  of  M  points,  the 
MSB  is  minimized  for  the  optimal  feedback  gain  .  The  optimization 

problem  then  appears  to  minimize  the  quadratic  performance  index  (cost  function) 


min 


“^{[^n+1  +  (^^|^_^))]  ]k<n<k+M-\ 

\<k<N 


]_ 

M 


k+M-\ 


N . 


(4a) 

(4b) 


Such  a  problem  may  be  solved  in  two  ways,  namely  by  employing  either  the  state  space  (HR)  model 
( A:  =  1 )  or  the  FIR  structure  {k  =  N).  Below  we  consider  both  cases. 
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STATE  SPACE  PREDICTIVE  ESTIMATION  OF  TIME  ERROR 


In  the  state  space,  the  clock  state  equation  and  the  observation  equation  are  given  by,  respectively. 


+  u„  +  ,  (a)  -  HX,„  +  ,  (b), 


(5) 


where  =  {Mintin' "^mn\  vector  of  a  control  signal  of  mxl  dimensions  (for  the  clock  model,  m  =  2  or 
3)  and  dimensions  of  all  other  vectors  and  matrixes  depend  on  the  number  of  the  states.  The  noises 
and  are  jointly  uncorrelated,  so  that  for  all  k  and  n  there  is  =  0 .  The  GPS  discrete  station¬ 
ary  white  Gauss  noise  has  zero  mean,  ^  known  covariance 


Ein^n^u) 


\V=a, 


On’ 


0, 


k  =  n 
k  ^  n 


(6) 


which  is  usually  assumed  to  be  time-invariant,  (7^^  =  V  =  (7^  =  N^/ 2 A ,  where  Nq/2  is  a  double-sided 
psd  of  the  corresponding  continues  noise.  In  the  clock  control  problem,  we  ignore  the  fact  that  the  noise 
is  colored  and  assume  it  to  be  white^.  Accordingly,  we  go  to  the  idealized  covariance  matrix 


(7) 


in  which  dimensions  of  T  depend  on  the  number  of  the  states.  For  the  sake  of  definiteness,  we  finish 
now  with  a  consideration  of  two  aforementioned  cases  of  the  state  space  filter  structures. 


Case  I — Two-state  filter  structure.  When  an  oscillator  noisy  phase  is  assumed  to  be  Brownian  with  a 
linear  drift,  only  two  states  become  non-zero  and,  thus,  the  state  transition  matrix  of  2x2  dimensions, 
measurement  matrix  of  2x1  dimensions,  and  error  matrix  of  2x2  dimensions  are  formed  by,  respectively. 


A  = 

1  A 

,(a),  H  =  [l  0],(b), 

"0  o' 

_0  1_ 

2 

_0  1_ 

(c), 


(8) 


where  N^/l  is  a  double-sided  uniform  psd  of  the  second  state  noise,  which  is  dependent  on  a  sample 
time  A.  We  notice  that  the  quantity  represents  the  desired  (or  conventionally  assumed)  signal- 

to-noise  ratio  (SNR)  in  the  observation  . 


Case  II — Three-state  filter  structure.  Assuming  the  oscillator  phase  to  be  Brownian  with  a  quadratic 
drift  leads  to  the  3-state  situation,  in  which  the  state  transition  matrix  of  3x3  dimensions,  measurement 
matrix  of  3x1  dimensions,  and  error  matrix  of  3x3  dimensions  become,  respectively. 


^  In  such  a  model,  the  noise  is  not  oscillator’s  white  phase  (f^)  but  rather  a  permitted  time  error  for  the  locked  clock 
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A  = 

’l  A  AV2" 

0  1  A 

,(a),  H  =  [l  0 

'0  0  O' 

0  0  0 

0  0  1 

2 

0  0  1 

(9) 


where  Nj^/ 2  is  a  double-sided  uniform  psd  of  the  third  state  noise,  which  depends  on  A.  Like  the  case 
(8),  here  the  quantity  represents  the  desired  SNR  in  the  observation  . 


The  linear  Kalman  estimation  algorithm  associated  with  the  filtering  and  1-step  prediction  problem  reads 


Enter  and 

Prediction  (project  ahead): 

(10) 

Compute  Kalman  gain: 

(11) 

Update  estimate  with  measurement  : 

(12) 

Compute  error  covariance  for  updated  estimate: 

(13) 

Make  1-step  ahead,  n  =  n+  1,  and  go  to  (10). 
Compute  1-step  prediction: 

(14) 

Predict  1-step  ahead  error  covariance: 

(15) 


R.=AR..,A"  +  T 


K.=R.H'(HR.H'+V)- 


1„  =  a1„.,  +  k.(4.-u.-ha1..,), 

R.=(I-K„H)R., 


^«+l|n  “  ^'^n\n  ' 


n\n 


where  in  (14)  is  the  estimate  (12)  taken  for  the  observed  n  points  and  in  (15)  is  the  error  co- 
variance  associated  with  the  estimate  . 


PREDICTIVE  FIR  FILTER  FOR  TIME  ERROR  ESTIMATION 

For  FIR  structures,  the  clock  time  error  equations  (1)  and  (2)  modify  to 

3-states:  ^  ’  (^)’  T.  =  To  +  ’  (1^)’  (c)’ 

(16) 

2-states:  ’  (^)’  Tn  =  To  ’  (1^)’  (1^) 


where  Xq  ,  yo  ,  and  D  are  initial  constant  values.  The  time  error  model  and  its  observation  equation  are 
then 


-  A„Xo  +  u„  +  ,  (a),  4  -  H>.„  +  ,  (b), 


(18) 
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where  for  the  constant  initial  state  vectors  and  kg  =  the  transition  matrixes  of  3x3 

and  2x2  dimensions  become  nonstationary,  respectively, 


’l 

A« 

0 

1 

An 

and  = 

1  An 

0  1 

0 

0 

1 

In  view  of  (18b)  and  =  x^  ,  the  FIR  filtering  estimate  of  a  time  error  reads 


N-l 


i=0 


z^W  =x^W 


(19) 


(20) 


where  Wi  is  a  weighting  function  (kernel  or  impulse  response)  of  a  FIR  filter;  z^(w)  =  is 

the  measurement  data  vector,  of  dimensions  Nxl;  =  [Wq(N),W^(N),...,Wj^_^(N)]  is  the  filter  weight 
matrix,  of  dimensions  A^xl;  xJ^{n)  =  {x^x^_^...x^_^^{\  is  the  time  error  data  vector,  of  dimensions  Nxl,  and 
\v^^{n)  =  [n^n^_^...n^_^^{\  is  the  GPS  noise  data  vector,  of  dimensions  Nxl.  Inherently,  the  first  estimate 
provided  by  (20)  appears  at  the  point  n  with  a  delay  on  the  transient 

^  =  A(A^-1).  (21) 


So  long  as  in  the  control  strategy  a  predictor  is  removed  out  of  the  FIR  filter,  a  1-step  linear  prediction  is 
formed  for  (20)  in  the  recursive  form 


^n-\-\\n 


^n\n  ^^^n-k\n-k  ^n-k-l\n-k-l^  ■ 


(22) 


where  an  average  increment  of  the  time  error  estimates  is  taken  for  the  reasonable  number  k  of  measure¬ 
ments  in  the  nearest  past.  Let  us  notice  that  there  may  be  utilized  some  other  predictive  algorithms  based, 
for  example,  on  the  regression  approximation  [1]  and  polynomial  FIR  predictor  structure  [6].  The  mean 
square  error  (MSB)  of  the  FIR  estimate  8^^  =  x^-x^  over  the  observed  database  is  calculated  by 

{£^  )  —  =  Nx^  +  ,  where  the  estimate  bias  and  variance  are  computed  by,  respectively. 


=  E{x„  -  ,  (a),  <  =  £[(x„  -  x;;W^  -  n;;W^  -  )^] ,  (b).  (23) 

We  notice  that,  in  contrast  to  the  state  space  estimator  whose  structure  is  exhaustedly  determined  by  the 
number  of  the  states  and  properly  filled  matrixes,  the  waiting  function  of  the  FIR  filter  is  a  matter  for  op¬ 
timization.  Below  we  consider  two  cases  suitable  for  the  GPS-based  locking  of  crystal  clocks. 

Case  III — FIR  structure  with  a  constant  kernel.  The  CK  FIR  filter  is  also  known  as  a  simple  moving 
average  (MA)  filter  that  is  optimal  in  a  sense  of  minimum  produced  noise.  Its  weight  is  rectangular,  i.e. 

where  =[1  1  ...1]^  is  a  unit  matrix  of  dimensions  Nxl.  A  disadvantage  is  a  large 

bias  when  a  time  error  function  demonstrates  an  easily  seen  nonstationarity.  Errors  associated  with  the 
filter  application  to  the  Brownian  phase  with  a  linear  drift  are  examined  in  [4]. 
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Case  IV — FIR  structure  with  a  linear  kernel.  The  LK  FIR  filter  is  optimally  unbiased  when  a  time  error 
function  behaves  as  a  Brownian  phase  with  a  linear  drift.  The  filter  was  designed  in  [7]  with  its  weight¬ 
ing  function 


2{2N-\)-6i 
N{N  +  \)  ’ 

0, 


0<z  <  A^-1 
otherwise 


(24) 


Statistical  errors  of  this  filter  for  the  Brownian  phase  with  a  linear  drift  are  discussed  in  [4].  An  advan¬ 
tage  is  that,  having  a  relatively  simple  kernel  (24),  the  filter  demonstrates  an  intermediate  performance 
between  the  three-state  and  two-state  Kalman  filters  and  this  is  independent  of  the  time  error  function. 


NUMERICAL  STUDIES  OF  OPTIMAL  LOOP  GAIN 

We  continue  on  with  numerical  studies  of  the  loop  (Figure  1)  employing  four  above-described  estimators. 
To  compare  the  results,  we  exploit  the  same  measurement  of  the  time  error  of  a  crystal  clock.  We  use  the 
Motorola  GPS  UT+  Oncore  Timing  receiver;  the  low-cost  crystal  oscillator,  AT-cut,  5  MHz  in  the  local 
clock;  and  the  time  interval  meter  with  a  resolution  of  1  ns.  A  sample  time  is  set  to  be  100s.  In  parallel, 
we  measure  the  time  error  for  the  reference  rubidium  standard.  We  numerically  solve  the  control  problem 
(4b)  for  the  steered  time  error  (1 -state),  assuming  the  oscillator  to  be  uncontrollable  (indirect  steering).  A 
typical  result  is  shown  in  Figure  2,  where  is  a  GPS-based  measurement;  \  is  an  actual  time  error 

measured  for  the  rubidium  standard;  is  a  controlled  measurement;  and  \  is  a  time  error  of 

a  locked  clock.  To  evaluate  the  trade-off,  we  set  the  same  time  constant  for  all  filters.  In  doing  so  for  the 
Kalman  filters,  which  has  a  IIF  structure,  we  set  the  noise  psd  Nyl2  in  (8c)  or  /  2  in  (9c)  to  finish 
the  transient  of  the  solution  of  the  discrete  time  Riccati  equation  (10)  at  the  level  of  0.95  with  respect  to 
the  first  state.  We  then  solve  the  control  problem  (4b)  for  the  optimal  feedback  gain  . 

Figure  3  shows  the  RMSE  of  the  locked  clock  as  function  of  k^^  for  several  numbers  N  of  the  points  in 
the  average,  employing  the  2D  Kalman  (a)  and  3D  Kalman  (b).  Figure  4  makes  it  for  the  CK  FIR  (a)  and 
LK  FIR  (b).  In  these  figures,  the  minimum  value  of  k^^  is  optimal  in  a  sense  of  LMS.  For  the  sake  of 
processing  accuracy,  we  refrain  from  setting  N  >  150,  having  a  limited  database.  We  notice  that,  by 
increasing  k^^  ,  the  clock  loop  evolves  from  a  large  error  to  its  minimum,  and  then  goes  through  the  in- 
termittence  to  divergence  (full  instability).  The  generalization  is  shown  in  Figure  5.  Here,  the  first  plot  (a) 
shows  how  the  optimal  gain  k^^  changes  if  N  increases,  and  the  second  one  (b)  demonstrates  the  same  for 

the  RMSE.  It  is  neatly  seen  that  the  unbiased  filters  (2D  and  3D  Kalman,  and  LK  FIR)  behave  with  quite 
some  similarity,  whereas  the  CK  FIR,  with  its  large  bias,  demonstrates  an  opposite  sign  of  the  slopes. 


DISCUSSION 

GPS-based  locking  of  a  local  clock  assumes  estimation  of  the  time  error  with  its  following  steering  in  the 
closed  loop.  Very  often,  the  estimation  is  provided  by  employing  FIR  and  state  space  structures.  To  de¬ 
termine  the  trade-off  between  these  filters,  we  simulated  and  studied  a  digital  control  loop  of  a  local  clock 
for  four  most  simple  filters:  2D  Kalman,  3D  Kalman,  CK  FIR,  and  LK  FIR.  Setting  the  same  time  con¬ 
stant  for  all  filters,  exploiting  the  only  typical  measurement  obtained  for  the  crystal  clock,  and  assuming 
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the  oscillator  of  the  clock  to  be  uncontrollable  (direct  discrete  control  of  a  time  scale),  we  arrive  at  the 
following  conclusions: 

•  With  small  TV,  when  the  noise  dominates  in  the  averaging,  the  CK  FIR  produces  better  result  over  all 
other  examined  estimators; 

•  With  large  N,  the  2D  Kalman  is  best  and  the  LK  FIR  occupies  an  intermediate  place  between  the  3D 
and  2D  Kalman; 

•  The  2D  Kalman  and  CK  FIR  have,  respectively,  less  and  more  sensitive  structures  to  variations  in  the 
optimal  feedback  gain. 

Overall,  we  notice  that,  for  the  sake  of  accuracy,  the  estimator  needs  to  be  designed  to  possess  a  structure 
that  guarantees  minimum  error  for  all  feasible  situations  with  the  clock  time  error.  Such  an  estimator  may 
be  produced  in  different  ways  [8].  However,  important  constrains  are:  1)  the  algorithm  must  be  simple  (it 
should  not  be  burdensome  for  clock  engineers),  2)  it  must  be  robust,  and  3)  it  must  rely  on  small  memory 
of  a  low-cost  crystal  clock. 
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Figure  1.  A  discrete  time  model  of  GPS-based  steering  of  a  local  clock  time  error  with  an  indirect 
control. 


Figure  2.  Typical  measurement  of  a  time  error  of  a  crystal  clock:  is  a  GPS-based  observation,  is 

an  actual  time  error  measured  for  the  rubidium  reference  clock,  means  an  observation  in  the 

control  loop,  and  x^  +  means  a  steered  time  error. 
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Figure  3.  RMSE  of  the  steered  crystal  clock  as  function  of  the  feedback  gain  for  several  numbers  N  of  the 
points  in  the  filter  transient:  a)  2D  Kalman,  b)  3D  Kalman. 
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Figure  4.  RMSE  of  the  steered  crystal  clock  as  function  of  the  feedback  gain  for  several  numbers  N  of  the 
points  in  the  filter  kernel:  a)  CK  FIR,  b)  LK  FIR. 
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Figure  5.  Statistical  performance  of  the  steered  crystal  clock  for  four  filters  as  function  of  a  number  N  of 
the  points  in  the  average:  a)  Optimal  feedback  coefficient  ,  b)  RMSE. 
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